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meridional wind component, positive south to north in m/ s. 


The components 0 and W define the polar form, and the U-V components define 
the Cartesian forms; 



Figure 1. The meteorological coordinate system. 



It Is more convenient to work in the Cartesian coordinate form for the 
subsequent statistical analysis than in the polar form. The statistical tools are 
readily available for use in Cartesian form, whereas for the polar form the 
statistical treatment would be awkward and in some cases techniques would havp 
to be developed. 
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Section IV, ’’Vector Wind Profile Models,” is the primary topic of 
interest in this report because of the potential applications to the flight per- 
formance analysis of aerospace vehicles. Several wind models are presented in 
this section. These models are called ’’synthetic vector wind profile models.” 
The concept follows that of the synthetic scalar (wind speed) model derived 
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distribution. The square root of the variance is called the standard deviation; 
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To emphasize the connotation of probability, F(t) is shown in Table 1 as P{x} 
versus selected values of t. The t-values in Table 1 are used as multiplier 
factors for the standard deviation to express the probability that a normally dis- 
tributed variable, X, is less than or equal to a given value. 



would be called the inter-95th percentile range. In common terms, 95 percent 
of the X’ s fall within Xi and Xo. 



TABLE 1. VALUES OF t FOR STANDARDIZED NORMAL (UNIVARIATE) 
DISTRIBUTION FOR PERCENTILES AND INTERPERCENTILE RANGES 



OJ' i'OOii' QUALITY! 














4 




b 


X 

1 

<M 



^ 0 
U 

o 'J-* 
^ CQ 
CQ O 

<D O 

s ^ 

0 X 


<0 Sj 

•B 

o 


<D O 

'S ^ 

al 

CD ?H 

i-rt ^-> 

^ W 

< 1 > *r^ 

1 -H 

m cd 

^ .s 

N bJ3 

u ^ 

S 5 

■s s 

u ^ 

oi is 

' « c 3 


•§ I § 

g ^ *43 

g a •? 

500 

g W "O 
cd *iH u 

^ fl ^ 
00 ^ 


+■> *r* 
OT bfc 


• 1-4 ^ 

^ -§ 

Sh 

cd 

o § 

Sh M 


^ Cd 
O T* 

^ b 


|Z! IH 0 

a 

rt § a 

a f-i TO 

uo 
o ^ ^ 

!=l O P 

^ g ^ 

o a 

kT ffi S 

O 

^ fl 

'S ® 

§ o "S 

^ ^ 

tM , , :^ 

M tH 

§ r I 

'-C 3 a ^ 

s a o 

d o +j 
+3 S 3 ra 

.2 ® 

^ 

I — I -i-i cd 
cd o 

•Sc .2 

Sh ^ -{3 

% ^ n 

a 42 m 


TO 9 h 


.2 2 

42 TO a 

^ ^ o 

■-S 1 

TO S O 

^ o 

cd ^ p 

s o o* 

r ^ 

^ cd 

O *E^ 

o > 

tJ *?3 !=! 

.2 a o 

s *43 

cd <D cd 

> a 

-M rr 
CD P qn 

^ CD 

^ TO <Z! 

p cd o 

^ o S 3 


^ p g 

-4-i I— 4 -4^ 

^ § I 

^ ” hc 

s ,a 

2 g ^ 

csS o (V) 


00 



o 



CQ 
• O) 

. CQ 

^ a 


m 

o X 

S 

o .5 
Pm o 
a 


l--^ 

^ -Sc 

1 -§ 

2< ^ 

o ^ 
o 


O fc£ 

-5 ^ 

m <1^ 

^ r/% 



where R(x) defines the region bounded by the ellipse of equation (9). The 
integration of equation (lO) is obtained by changing the variables X and Y to 
polar coordinates. The result is 


Solving for 7^ and replacing P(x) by p, we have 




TABLE 2. VALUES OF X FOR BIVARIATE NORMAL DISTRIBUTION 

ELLIPSES AND CIRCLES 





Because graphical displays of the probability ellipses are useful and 
informative, we have devised a convenient plotting procedure for electronic 
computer operations. First, rewrite equation (9) in the more recognizable 





For graphical presentations the range of the variable is important in 
order to arrange the scale. The largest and smallest values of X and Y for a 
given probability ellipse, p, are given by 
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ellipse. From this point of view, a specified probability ellipse gives the joint 
probability that p-percent of the U-V components lie within the given ellipse. 




Values for X for selected probabilities, p, are given in Table 2. 


Because this function is simple, it can be easily graphed manually. 
However, the generalized plotting technique for electronic computer plotters as 
represented by equation (14) can be advantageously used. 



(univariate) distributed. The resolution of bivariate normally distributed 
variables (the modulus of the vector) is Rayleigh distributed. The sum and 



differences of bivariate normally distributed variables are normally (univariate) 
distributed. For bivariate normally distributed variables, the probability distri- 
bution remains bivariate normal and the marginal distributions remain normal 
(univariate) with the rotation of the coordinates. 
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The asterisk is placed on X (i. e. , X*) to emphasize that X is the given value 
which may be assigned any arbitrary value. For illustrations and many applica- 
tions, it is convenient to let X’*' = X + to , where t is chosen from Table 1 to 
give a desired probability level for X*. 




and conditional standard deviation 


The above conditional normal probability distribution functions are univariate 
normal distributions for a ( fixed) given value for one of the bivariate normal 
variables. Thus the t-values given in Table 1 are applicable for conditional 
probabilities statements. For example, 
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This conditional distribution is bivariate normal. After computing the two con- 
ditional means, the two conditional variances, and the conditional (partial) 
correlation coefficient from the required 14 parameters describing the quadra- 
variate normal distribution, the joint conditional probability distribution function 



is computed using these conditional parameters as inputs to the same equations 
used for the bivariate normal probability distribution function of Section II. B, 
Although the extended algebraic expressions, even for this simple case, for the 
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the major and minor axes of the bivariate normal probability ellipse. 
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and other special cases are discussed by Smith [4] and more recently by White 
[7] , who gives moments for some special Rayleigh functions. The density 
function, equation (29) , can be integrated in closed form over any range of the 
variable R. Hence, the Rayleigh probability distribution function, F(R), is 


b. Bivariate Rayleigh Distribution. The mathematics for tlie gen- 
eralized bivariate Rayleigh probability density function become very complex. 
Following the work of Miller et al. [8] for the special case where the variables 
Xj and X 2 are bivariate normally distributed and independent with a = cr = 
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defined. For the special case when x = y = 0, equation (33) reduces to the 
following simple case: 
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The probability density function for the difference of two variables rj 
(X-Y) for the stated conditions is normally distributed with 
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This is the usual rotation of variances found in texts on statistics and probabil- 
ity. X (T and X (T , give the major and minor axes for the probability ellipses 
e a e b 

Fx = '^2 ^ -In (1 - p)l. 
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Rotation of the linear correlation coefficient through a degrees; 



where cov (X, Y) is the rotated covariance, 
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Note this change of order for the sign of U’ and V’ as H cross H . This is 
important later in computing the vector wind profile model. 
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we have results for the generalized Rayleigh distribution computed by equation 
(28b), where the sample parameters are u, v, s , s , andr(u,v). An example 




for the March wind speed over Cape Kennedy at 4, 8, and 12 km altitude com- 
pared with empirical percentiles is shown in Figure 7. Many other examples 
giving the results of percentiles for wind speeds derived from the generalized 
Rayleigh probability distribution all agree very closely with the corresponding 
empirical percentiles. 
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TABLE 3. FREQUENCY OF WIND DIRECTION AT 12 km ALTITUDE, 
FEBRUARY, CAPE KENNEDY, FLORIDA 
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The period of record for, g(o) , and observed fre- 
frequency is from 1956 through 1967. There is 
excellent agreement between the derived frequency 
of wind direction and the observed frequency in this 
example. 




TABLE 4. FREQUENCY OF WIND DIRECTION AT 10 km ALTITUDE 
DECEMBER, VANDENBERG AIR FORCE BASE, CALIFORNIA 
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0.4026. The observed frequency is from the data sample 
from 1965 through 1969. Considering this, the agreemeni 
between the derived frequencies of wind direction and the 
observed frequencies is excellent. 





I s I 

^ 'g o 
a ^ a 


(D 

'TJ be 
Cd a 

S3 

^ S 

o ^ 

g J 

^ cd 

a ^ 

.2 CQ 

■e 0) 

^ s 

*»^ ^ 

;s 
M Id 
-3 ^ 


be ^ 

Ch 'S 

•rt .2 

M cd 
3 X! 

m ° 

^ a 

•S rt 

TS O 
T3 

S 

^ § 

o o 

O Q 


cd 

•s © 

S-I 

s 

” s 

^ d) 

w .S 

• cd 

1 — I ►> 


cd 


a 

a; 

a, 

a 

o 

o 


cd 
o 

Q 

O 
O 

cd Q 

be 

c o 

•rH r— t 

cQ a^ 

o ^ 
O OQ 

?n Cd 

be ^ 

cd cd 

s ^ 

cd <D 

5 

OT 

•+-> “S 

c^ a 

0 o 

be o 


0 


cd 

a S 
© ^ 

Eh (X) 

cS 

© 

N © 

a 

o S 


0 

>H 4-» 

p 

u ^ 


p a 0 

Eb p 


0 

H 


0 

cd 

& 


o ^ 
0 w 

-H 

CO B 
CO ^ 

a 2 

O rC 

•43 

s a 

^ I 

.2 CO 

^ "S 
0 

a 

o 

a 


0 

> 

0 

4^ 

o 


rC 

0 

I 

0 

~Su 


cd 
0 

•a 

u 

a 2 

0 ^ 


o 

0 


0 

r-H 

Oh 

a 

cd 
CO 

0 
XJ 

sd 
0 
> 

•iH 

be 

§5 

§ 

u 

T3 
0 

t 
0 
CO 

-8 

0 
T3 

•iH 

CO 
-JJ 

§ 
T3 
0 

•iH 

^ -3 

P o 


cd 
CO 

CO 

cd 2 


CO 

0 

Id 

a 

4J 

CO 

0 


•8 

-8 

U 

Oh 

CO 

§ 


-o ■« 

I ^ 
I 

S w 


"0 


lO 


Q 

W 

W 

^ 

c/2 

Q " 
2: ® 


o 

>1 

Eh 

I— I 

>-i 

(— I 

pq 

< 

m 

o 

A 

hPi 

< 

o 

HH 

Eh 

I— I 

Q 

2 

O 

O 

lO 

H 

hJ 

pq 

<3 


2 

O 

HH 

H 

U 

w 

Q 

Q 

:z; 


< 

iz; 

w 

s 

o 

CO 


lO 

i> 

VI 

* 

o 


g 


w 

:zj 

:z; 


w 




H 

pi< 

< 

u 



(N 

tH 

H 

< 






225 SW 76.79 



^ 05 

CO CO 


LO <M 
OC C75 


O LO 
l> 'tH 
eg CO 



ZD 




CO 

L- 

cu 

CD 

-C 

oo 


o 

CD 


"O 

c: 


o 

CO 

c: 

o 

H— ' 

13 

-O 

"Co 

s 


o 


to 

c3 

0 

m 


?-< 

O 

4-i 

o 

0 

> 

o 

§ 

+-> 

U1 

;g 

3 

'u 

> 


0 

H 


•bJD 

CQ 

Pi 


O 


O 

§ O 
o > 

H •$ 

<N ^ 
pi 
O 

+3 

u 

C3 

a 


0 

B 

-o 2 

PI cJ 


u 

O 

4-» 

o 

> 

«4-l 

O 


0 

G iS 
O ^ 




•§ 

U 

Oh 

rG 

CO 

G 


O 

0 . 

> O 

«4H 

o 

G 0 

0 T3 

^ B 

1 S 

U cj 

*CQ ^ 
.G 0 

0 

U 

S 

B ^ 
^ :S 

i 'S 

^ «) 

G 

O 


G 


0 

I-H 

a, 


G g 

.2 ^ 
0 

0 CO 

■s S 

•fH *+^ 

^ O 

> O 

2 ^ 

0 . 

B 

CO 2 

2 I 
0 0 
«5 lijd 
0 

G 0 
a Oh 

0 ^ 


I 

X 

0 

CO 

G 


u 

G 

> 

O 

o 

G 


0 

H 


G 

a 

CO 

G 

CCS 

0 

rG 

CO 


G 

O 

•4-> 

0 
0 
> 
•+^ 
G 
0 
0 
G 
0 
G. 

05 

05 

§ 

4-> 

0 

0 

G 

0 

Gh 

LO 

05 

0 

T3 

§ 

1 

(M 


CO 


^ ffi 

M 

G 0 

HH ^ 
0 § 

•2 13 

'rt t 

^ 0 

2 2 


a 


2 ^ 
H G 

0 
g: 

E 

0 

r-l 0 

> 

ci *bO 

•I— I 

.2 

— ^ G2 

^ H 

§ s 

CO G 
G 05 

si 

«« o 

> *413 

» t— H 

^ G 

Ph 0 

02 

CO 

0 0 
CO 

a 

Ph g 
0 il 


s 


CO 


0 

G 


5 


0 

G 

CCS 


0 

> 

o 


G 

> 

G 

0 


0 

G 

a 

•tH 

Pg 


CCl 

CO 


T3 


jy 2 

_ .rG 

S H 

B 


CO 


ffi 

I 

c 

m 

§ 


1 

CO 


<M 

I 

CO 


00 

I 


G 

G 

O 

h-* 

HH 

hIh 

HH 

0 

X 

CO 

T3 

G 

1 

T3 

G 

G 

HH 

G 

2 


P. 


- G G G 

00 

X 

S 

X 

T3 

00 


LO 

G 

tH 

tH 

tH 

G 

CO 

CO 

CO 


•fH 

•fH 

•iH 

(M 

K 

ffi 

ffi 


I I 


T3 

0 


tH 

' T— 1 

rH 

T— [ 

tH 

tG 

0 

0 

0 

G 

G 

G 

§ 

a 


s 

s 

E 

• 

• 

• 


G 

'TT 

G 

G 

•> 

co" 

co'' 

•« 

1 > 

»• 

I G 


CO 

G 

0 

•G> 

0 


G 

G 

G 

G. 

G 

G 

0 

.G 

CO 

T3 

G 


G 

O 

0 

0 

> 

> 


0 

X3 

H 


i T5 
G 

w g 

•rH ^ 

f ■ I 

•fH ^ 
fG O 

G CQ 

I « 

G. .& 

I— H 

G 0 

S 

O G 

.5 o 

G ^ 
G a 


(M 

lO 


I> 


G 

G 2 
G 

O w 
Gog 

rG 0 

^ Cu3 r5 

G 


0 

> 

0 

G 

G 


ro 

G 


CO 


S § 

0 X 


^ fS 

a a 

CO 

G g 

*43 o 

G G 
& 

O G 
r-( O 
G 0 

§2 

0 ^ 
bJD 

«c 

G 

^ .2 
•iH 43 

0 

CO 

^ 0 


a 

CO 

G 

O 


G 

O 


O 

0 ^ 

^ & 0 
C G 0 

o *'■" > 

G oS O 


^ s ^ 

0 2 G 

X G 
CO 1^ ^ 
TO G 

S 9 P 


G 

0 

G 

0 

X 


0 

G 

2 

W 

0 

4-> 

cj: 

G 

•rH- 

I 


G 

O 


CO 
0 
CO 

2 U 

-IS 

o o 
5 CO H 

O H 03 
^ G 

. 0 

G 

® I fl 

G o G 
^ ^ G 


(?j X CQ 0 


Brooks and Carruthers [13], was applied to these sample parameters: 
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Figure 11. Bivariate normal distribution of vector wind shears (m/ s) , 

F(u’ , v’ )» reference altitude H = 12 km and H = 9, 10, 11, 13, 14, 

o 

and 15 km, March, Cape Kennedy, Florida. 
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Figure 13. Conditional bivariate normal distribution of vector wind 
shears (m/ s) , given vector wind at = 10 is 330 degrees, 58 m/ s, 

F(u» ,v* |u* = 28.9, v*= -50.1), H= 8, 12, and 18 km, December, 
H H H 

*^Vandenbere Air Force Base, California. 
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by the desired scientific or engineering application. In Section IV. B a step-by- 
step procedure is given to compute the SVWP that is in-plane with the given wind 
vector. This in-plane profile has two branches: one is the smallest conditional 
vector wind and has the largest shears, and the other is the outer branch, which 
has the largest in-plane conditional wind vector but not necessarily the largest 
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Figure 15. Conditional bivariate normal vector winds given the wind 
vector at = 10 is 330 degrees at 58 m/ s, December, Vandenberg 
Air Force Base, California. 



Select the probability, P , for the probability ellipse. 
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cover all quadrants; however, the range of 0* that will intercept the ellipse 
is required. To determine whether the probability ellipse occurs in all quadrants 
for a given probability ellipse, compute 




Figure 16 provides definitions concerning the derivation. The requii 
ment is to obtain the lines that pass through the origin and are tangent to the 
ellipse, F(x,y), 



Figure 16. Definitions for line intercepts to an ellipse. 
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The wind vectors {Oj ,WiJ and {02 > W 2 / are the tangent intercepts to the 
probability ellipse at H . 




The slope, m, is computed in this way to preserve the meteorological sign 
convention. 
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Compute the expected values: 
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Note that the conditional standard deviations are a function only of the 
correlation coefficients and shear standard deviations, and these parameters 
are a function of and H. We are using the simplifying assumption that the 

conditional vector shears can be treated as bivariate circular distributed. 
Therefore, compute the resultant conditional standard deviation 
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E(v|v*) = V* + E(v’|v*) 



At H = H , E(u* |u*) = 0, hence E(u|u*) at H = u*. Similarly, at H = H , 
o o o 

E(v’ |v*) = 0, and E(v|v*) at = v*. The change of sign between equations 

(66) and (67) is because of the convention used in computing the wind component 
shear sample values [see equations (50) and (51)]. 
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both solutions with respect to the 0*-plane. Compute 


and take the smaller of the two values as W 4 . If Wi is from the interce 
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When Wj >0, 9\ = 0*; when Wj < 0, 0j = (0* - 180°) . The larger intercept, W 

A 

from equation (71) is also an in-plane SVWP but is of the same direction as 0*. 
This profile can, and often does, exceed the specified wind vector probability 
ellipse. 
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direction from which is blowing, i. e. , the meteorological convention. Also 
note the following: 
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Note that the quantity A. s is the radius of the conditional vector wind probabil- 
ity circle. ^ 


c. Condition C. When Wj < 0, and E(u|u*) and E(v|v*) have opposite 

signs to u* and v*, and (W - s ) <0, the conditional probability circle has 

E Vj[ 

passed to the opposite quadrant from u* and v* and is completely contained in 
the opposite quadrant. This condition will be rare except for P « 0. 95. Th( 
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The curve labeled W^ is a valid SVWP that is in -plane with 0*, but it 

exceeds the 95th percent wind ellipses at all altitudes except at the reference 
altitude H^. Some rationale to exclude this case based on an evaluation of the 

SVWP concept for aerospace vehicle response to wind should be established. 


The wind shears obtained from SVWP (2) and SVWP (3) (Figs. 18 and 

19) are less than those for SVWP (l) (Fig. 17) . The vector wind shears for 

SVWP (l) for this example {O* = 330 degrees, W* = 57.8 m/s, II = 10 km, 

o 

December, Vandenberg Air Force Base) are much less than those of Reference 
14. However, for other given wind vectors (o*, W*} for SVWP (l) , the in-plam 
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Figure 17. Synthetic vector wind profile, SVWP (l) , in -plane with given 
wind vector, December, Vandenberg Air Force Base, California. 
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Figure 19. Synthetic vector wind profile SVWP (2) , wind vectors tangent 
to the left of conditional wind circle, December, Vandenberg Air Force 

Base, California. 




It must be recognized that comparisons are being made between two different 
wind modeling concepts and that the statistical parameterizations are different* 
The wind characteristics over Cape Kennedy and Vandenberg Air Force Base 
are different; e. g. , the monthly vector mean winds over Cape Kennedy are 
greater than those over Vandenberg Air Force Base, and the vector wind 
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A method to obtain the SVWP with respect to tlie monthly mean vector 
wind is presented here. Another alternative is to obtain the conditional wind 
vector probability ellipses rather than the conditional circles. A summary 
guide for SVWP procedures completes this section. 
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When the vector wind probability ellipse is specified, e. g. , P = 0.95, 
the intercept of the conditional vector residual, {r*, <^*} to the probability 
ellipse can be obtained by simple calculations (Fig. 21) . Hence, with respect 
to the original coordinate system the wind vector that intercepts the 




probability ellipse is readily computed in terms of the components u*and v* 
This procedure assures that a given wind vector { 0 *, W*} at the reference 
altitude will intercept the probability ellipse. This method is less complex 
than that used in Section IV. B. There is also an advantage in this procedure 
in the application of the SVWP for aerospace vehicle trajectory control system 
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b. Alternative 3: SVWP with Respect to Vehicle Flight Azimuth. 
The SVWP may be determined by computing the intercept of the wind vector 
from the origin to the conditional probability circle where the direction is 




Vectors 1 and 4 are intercepts to the conditional vector wind 
circle in-plane with the monthly vector mean wind at the refer- 
ence altitude H , or versus II. 
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The tabulations of the 14 parameters for the quadravariate normal 
probability distribution of wind vectors and wind vector shears have certain 
symmetry for altitudes above and below the reference altitude H . This 

r\ 
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Space Division for suggestions leading to these alternative SVWP models. 




d. 99 Percent Vector Wind Shear. The 99 percent vector wind shear 
for the conditional probability of wind vectors given a wind vector on the 95 
percent probability ellipse should be used. 






h. Interpolation for Shear Interval Less Than 1 km. To obtain the 
SVWP over the the shear interval between reference altitude and 1 km below 
and 1 km above the reference altitude, use interpolation equation 5. 26 (page 
5 . 96 ) of Reference 14 . 
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Falls, L. W. ; Normal Probabilities for Vandenberg Air Force Base 
Wind Components — Monthly Reference Periods for All Flight Azimuths, 
0- to 70-km Altitudes. NASA TM X-64897, January 1975. 




Daniels, Glenn E. and Smith, Orvel E, : Scalar and Component Wind 
Correlations Between Altitude Levels for Cape Kennedy, Florida, an( 
Santa Monica, California. NASA TN D-3815, April 1968. 
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A PROGRAM TO COMPUTE CONDITIONAL BIVARIATE 
NORMAL PARAMETERS 
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the symmetric variance-covariance matrix. 



A property of the multivariate normal distribution is that marginal and 
conditional distributions are also normally distributed. The general expression 
for these distributions is found often in the literature. ^ Remarks are confined 
here to the specific case. 



Morrison, D. F. (1967): Multivariate Statistical Methods . Wiley, N. York. 
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The 5 parameters in the first row and the 9 parameters in the body of 
the table contain the 14 statistical parameters required for the quadravariate 
normal probability distribution function for the vector wind at a reference 
altitude, HO, and the vector wind shear at altitudes below and above the reference 
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